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Abstract 

We study the effects of Lorentz symmetry violation on the scalar CMBR bispectrum. We deal 
with dispersion relations modified by higher derivative terms in a Lorentz breaking effective action 
and solve the equations via approximation techniques, in particular the WKB method. We quantify 
the degree of approximation in the computation of the bispectrum and show how the absolute and 
relative errors can be made small at will, making the results robust. 

Our quantitative results show that there can be enhancements in the bispectrum for specific 
configurations in momentum space, when the modified dispersion relations violate the adiabatic 
condition for a short period of time in the early Universe. The kind of configurations that are 
enhanced and the pattern of oscillations in wavenumbers that generically appear in the bispectrum 
strictly depend on the form of the modified dispersion relation, and therefore on the pattern of 
Lorentz violation. These effects are found to be distinct from those that appear when modelling 
very high-energy (transplanckian) physics via modified boundary conditions (modified vacuum). 
In fact, under certain conditions, the enhancements are even stronger, and possibly open a door to 
the experimental study of Lorentz violation through these phenomena. 

After providing the general formulas for the bispectrum in the presence of Lorentz violation 
and modified dispersion relations, we also discuss briefly a specific example based on a healthy 
modification of the Corley-Jacobson dispersion relation with negative coefficient, and plot the 
shape of the bispectrum in that case. 
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1 Introduction 

The non-Gaussian features of the anisotropics of the Cosmic Microwave Background Radiation (CMBR) 
have received increasing attention in recent times, in connection with the upcoming release of the 
data of experiments like Planck [1]. The primordial scalar non-Gaussianity occurs when the three- 
and higher-point functions of the comoving curvature perturbation £ are non-zero, and are therefore 
a manifestation of the interactions at the time of inflation [2]. However, they could carry signatures 
of physics at much higher energy, in particular of the so called transplanckian physics 011]. 

In this work we are interested in the bispectrum, which is obtained from the three-point function, 
in single-field slow-roll models of inflation. In the conventional scenario this is slow-roll suppressed [2J, 
hence any non-Gaussian signature is a smoking-gun for deviation from the standard picture. The 
effect one would generally be interested in for detection is a peculiar shape function, that is a specific 
shape for the graph of the bispectrum as a function of the external momenta, with, in particular, the 
presence of enhancements for some of their configurations. 

There have been several approaches to the study of high energy (transplanckian) physics in cosmol- 
ogy: from the standard program of including higher-derivative terms in the field theory action, to the 
modeling via modified uncertainty relations [5], or via specific boundary conditions imposed at certain 
times/energy scales through boundary actions/momenta cutoffs j6j[7]. In this paper, we focus on the 
effects on the bispectrum due to higher-derivative corrections yielding modified dispersion relations for 
the perturbations fields 0- The dispersions differ from the Lorentzian one when the physical momenta 
are above a certain energy scale A 3> H, where H is the Hubble rate at inflation |1UH12| . 

Such effects would represent a signature of Lorentz violation at high energy and possibly open up 
new opportunities to study those phenomena through the CMBR. The interest in Lorentz violation 
has always been high [13H17] and revived in recent times also by certain realization of quantum gravity, 
such as the Hofava-Lifshitz model [18]. Lorentz violation also arises in braneworld models [19] and 

1 Our scenario differs from that considered in [8][9], where the oniy modification to the fieid equations was a change 
in the speed of sound (the modifications at the ievei of the Lagrangian depended only on first order derivatives of the 
fields). In particular, we include also derivatives of higher order, and the effect on the field equations is more profound. 
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modified dispersion relations also appear in effective theory of single field inflation when the scalar 
perturbations propagate with a small sound speed [20] . 

In |21| . the effects on the bispectrum due to modified dispersion relations which do not violate the 
WKB (adiabatic) condition at early times were studied via the specific example of the Corley- Jacobson 
dispersion relation with positive quartic correction to the momentum square term [22j, where the exact 
solution to the field equation could be obtained. No large enhancement factors were found, but the 
leading modifications to the standard slow-roll results were strongly suppressed by the ratio jp. This 
can be explained in terms of the very small particle creation, due to the absence of WKB violation at 
early times. 

In this article, we provide a more general analysis of modified dispersion relation and their effects 
on the bispectrum. In particular, we consider the case where the adiabatic condition is violated in the 
early Universe for a short period of time. In this more interesting scenario, particle production is more 
substantial. In fact, the Bogoljubov coefficients accounting for particle production do not depend on 
small ratios of energy scales such as but, as we will show, are only constrained by backreaction. On 
top of that, particle creation might also lead to enhancements factors in the bispectrum for specific 
configurations of the three external momenta, due to interference patterns. 

We will find that the leading modifications to the standard slow-roll result for the bispectrum 
in the presence of this kind of modified dispersion relations exhibit two main features: there appear 
modulations (oscillations) in function of the momenta and there can be enhancements factors for some 
of their configurations. We show how their presence, their magnitude and the enhanced configurations 
depend on, and can be obtained from, the form of the modified dispersion relation. We also find that, 
under certain conditions, the enhancements can be actually larger than those present within the mod- 
ified vacuum framework for transplanckian physics [3], in particular when there are higher-derivative 
interaction terms. The reason is that in the case of modified dispersion relations some particular 
cancellations in the three-point function do not occur any more. We comment on the likelihood of 
satisfying the necessary conditions for having enhancements. We also study the backreaction issue 
and derive the relevant constraints descending from it. 

The results are robust, as we will show that the errors deriving from the use, when necessary, of 
approximations to the solution of the field equation can be made small at will. Our general results 
can be easily adapted to the theory models of interest, obtaining from the effective action the relevant 
dispersion relation. In the appendix at the end of this paper, we apply our analysis to a particular 
example, which is a physical modification of the Corley-Jacobson relation with negative coefficient, 
avoiding imaginary frequencies. 

The outline of the paper is as follows: we start in section [2] by introducing the formalism for 
describing Lorentz violation and cosmological perturbation theory, and setting the notation. We then 
study the modified field equation for the perturbations and its solution in section [3j Our results 
regarding the bispectrum are presented in section [3J in particular, in section H~T1 we study the case of 
pure cubic interactions and in section 14.21 the case of cubic terms with higher derivatives. We then 
discuss the backreaction issue and find the constraints it imposes in section [5j We conclude in section [6] 
with a summary and comments. Finally, the typical Lorentz-breaking action can be read in appendix 
IA.H and we apply our results and techniques to a specific example of modified dispersion relation in 
section IA.21 

2 Formalism and notation. 

The approaches to the implementation of Lorentz violation have been various (for reviews, see |13j). 
We will prefer the most conservative point of view, within the framework of Effective Field Theory and 
preserving general covariance and absence of fixed geometry. These latter require that Lorentz sym- 
metry be violated by dynamical Lorentz violating tensors, which, if we preserve rotational symmetry 
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in the dispersion relation, can be reduced to (products of) vectors [13j . Preserving general covariance 
is clearly appealing, because otherwise we would have to give up general relativity. 

The standard procedure of cosmological perturbation theory can then be followed. It begins dis- 
tinguishing background values and perturbations for the various fields. The cosmological background 
has a scale of variation much lower than those of the Lorentz breaking corrections, and its standard 
description is valid. The distinction between background and perturbations can be operated using 
different threadings and slicings. In [2], two gauges were used, differing in the behaviour of the co- 
moving curvature perturbation t) and the inflaton perturbation ip{x, t) (we indicate with (j)(t) its 
background value). In one gauge 

C(x,t) = 0, <p(£,t)^0, (1) 

in the other 

C(x,t)^0, <p@,t) = 0. (2) 

The latter gauge is more convenient conceptually, as one works directly with the physically interesting 
perturbation (the curvature one, which is conserved outside the horizon), but the former one is often 
better computationally. 

In our case, the first gauge is by far the most convenient, as one can take into account the higher 
derivatives correction via the action of the inflaton (a simple scalar field). The inflaton sector will 
therefore be our Lorentz-breaking sector. This is especially appealing when discussing backreaction, 
as we will do in section [5l because one does not have to compute the effective stress-energy tensor of 
the metric perturbations, which is a complicate calculation. A covariant action with higher-derivative 
terms and a dynamical Lorentz-breaking vector field can be constructed [12 ^ 14 4 [T5]: we report it in 
appendix lA.li 

The transformation between the two gauges is a time-reparametrization t — > t + T, such thatH 
( = HT+ ±H >t T 2 - ^dtTdjTa" 2 + X -d iX diT + -^d^d^ (diTdjT) - ^cT 2 ^ {dadjT) (3) 
where H and a are, respectively, the Hubble rate and the metric scale factor, while 

rr _ V 1(^ 2 + W d\ = ^£U-^\ (4) 



</>, t 2 (<^) 3 (<A, t ) 2 ' m dt x Vit j 

12 

The higher order terms in ip and in the slow-roll parameter e = 2Hi ^ are necessary to compensate 

Planck 

for the time evolution of ip on superhorizon scales and make £ conserved outside the horizon. For 
our needs, it suffices to consider this relation only at leading order in perturbations and slow roll 
parameters. 

We expand the inflaton perturbation (p as 

/d 3 k 
—-j n (ri)e lk - x , (5) 

using the conformal time r) = f dta _1 (t), and quantize writing 

^ k (v)=f R (v)&l+fim % , (6) 

where fAjj), fHv) are two linearly independent solutions of the field equation 



fl+{"(v,k?- Z -)f- k = z= a -f. (7) 



^Derivatives with respect to t are indicated with the notation ,t, those with respect to the conformal time r\ with '. 
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Different choices for /|(??),/^(^) correspond to different choices of vacuum for the field. Imposing 
the standard commutation relation on the operators cL f , tL , implies a certain normalization for the 
Wronskian of f?(v)) fki 7 ])- Using equation Q at leading order, the comoving curvature perturbation 
is then expanded as 

J (27TJ2 Z 

In equation (J7|), w(ry, k) is the comoving frequency as read from the effective action. In the standard 
Lorentzian case w(ry, k) is equal to k = \k\, but for a modified dispersion relation oj will have a different 
dependence on k, rj. For isotropic backgrounds, which we will limit ourselves to, fu(v) an< ^ 
depend only on k and, therefore, we may drop the arrow symbol in the following. 



3 Modified dispersion relations, field equations and WKB method. 

We rewrite here for convenience the field equation ((7]) for the mode functions: 

fl + V(ri, kff % = V( V , kf = u>(tj, k) 2 - -. (9) 



Solving it in the case of a modified dispersion relation is often difficult and approximation methods 
have to be employed. One of them consists in exploiting the WKB approximation where possible. This 
approach provides a global approximated solution to differential equations whose highest derivative 
term is multiplied by a small parameter that we call e [23]. We will therefore rewrite equation ([9]) in 
such a way that the WKB method can be rigorously applied. Our approach will somehow differ from 
what has been done in the past in the cosmological literature discussing the spectrum of perturbations, 
and will enable us to have a better control on the approximation. 

It is useful to start with a brief general review of the WKB method to make the paper self-contained 
and fix the notation; more details can be found in [23]. The method applies to equations of the form 

d n f d n ~ x f d n ~ 2 f 
e-r^ + a(y)-, J T + b{y)-, £ + • • • + m(y)f = , (10) 

but it can also be used for inhomogeneous equations, to approximate the Green functions. In our case 
the equations are of second order and can be put in the form 



dy 2 



e 2 -± + q(y)f = (11) 



The WKB method consists then in: 

i) postulating an approximated solution of the form 

/(y)~e* (12) 

for a small parameter 5 that will be fixed by studying the differential equation as we will show, 

ii) substituting the ansatz (fT2|) in the equation (fTT|) . obtaining 

2 2 2 

^sl + y2s si + js + e 2 s 2 + e 2 2s s 2 + e 2 h + 0{e 2 5) + q(y) = , (13) 



3 From now on we indicate with dots the derivatives with respect to the variable y for which the differential equation 
at hands has the form 1)1111 . The apex ' will keep denoting the derivation with respect to 77. 
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Hi) determining the value of 5 through a distinguishing limit (see [23] and in the following) and 
solving the sequence of equations obtained from (|13p by comparing powers of e. 

We will soon give practical examples of these steps when presenting the cosmological application. Be- 
fore that, let us discuss the validity conditions for the WKB approximation. The first condition is that 
the series at the exponent of (|12|) is an asymptotic series uniformly for all the interval of x considered. 
Moreover, the validity of the procedure for solving the differential equations requires/guarantees that 
the series can be differentiated term-wise (at least up to the order of the equations) and the derivative 
series are all asymptotic as well. We then need to satisfy 



5s n+ i ~ o(s n ) 



5s 



n+l 



<C 1, , Ss n+1 ~ o(s n ) ->■ 



5s 



n+l 



« 1, 



Vn, 



(14) 



uniformly for the interval of x considered. 

Furthermore, since the asymptotic series is at the exponential, these conditions are not enough to 
guarantee that we have a good approximation of f{y): for the WKB series truncated at the order m 
to be a good approximation of f(y) we must also require 

5 m s m+1 < 1 5^0. (15) 

The relative error made when using the WKB approximation truncated at the order m is then 



f(y) - e * ^™=° 5ns " 



f(y) 



5 m s 



m+l 



(16) 



We will also be concerned with the error done when integrating the WKB approximation in place of 
the exact solution. It can be found that 



1 v-^ m A" 



Sm+1 



J J i=l J i=l j=l 



m+l 



(17) 
(18) 



1 v-^m cn 

which, beside being small because proportional to 5 m , is further suppressed when e*^ n =° i,n is a 
rapidly oscillating function, as it will be in our case. These criteria are quantitative and tell us the 
order we need to go to for approximating the result to some prescribed error, which we can make 
small at will by going to higher orders in the approximation. 

The conditions (|14p . (|15p are the WKB condition(s). Given an equation of the form (|lip . they 
generically require a slow variation of the quantity q(x) and are also called "adiabatic". 

Let us now investigate the application of the WKB method to the cosmological equation ([9]) in the 
case of modified dispersion relations. We will show that there is in fact a natural choice of variables 
such that a small parameter appears. We consider dispersion relations where the usual linear behaviour 
receives non-negligible corrections when the momentum is larger than a certain physical scale A 3> H, 
and therefore the dispersion relation can be written in general as |f| 



phys 



(P) 



P 



a(rj) a(r?) 



F(x -> 0) -> 1 . 



(19) 



4 In many cases, the dispersion relation is given directly as an expansion in power series, for example u)(k,ri) 2 » 
5Zn=o c "{%r) n (k 2 )" +1 r/ 2 " , Co = 1; we however prefer to use the generic formula in (|19[) . There are some discussions 
about the scale suppression of the higher derivative terms, see |13| . 
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We now change variable to y = kr/j^ and calo = e. The field equation then reads 

e 2 f y +(F(y) 2 -e 2 ^\f y = 0, (20) 

So far, we need not specify if we discuss sub- or superhorizon scales. We proceed with the WKB 
method and write f)13j) with 

q (y) = F{y) 2 -e 2 \. (21) 

yZ 

By looking at equations (|13p and (|21|) . we spot two possible distinguished limits for 5, that is 5 ~ e and 
<5 ~ 1. They arise from two different dominant balance conditions [23]. In the first case, the leading 
equation when considering powers of e is 

s 2 + F(y) 2 = (22) 
In the second case the leading equations are 

^ n n y ' * 

Here, it surfaces the distinction between superHubble and subHubble modes, which corresponds re- 
spectively to the first and second distinguished limit for S. In section I3.1| we will use the WKB 
approximation in the intervals of time where the first limit is appropriate, and look for a different 
solution in the case of the second limit. In fact, in this latter case, the WKB conditions (|14p . (|15p will 
not be satisfied H. 

For the moment, let us discuss which kind of modified dispersion relations can appear. In general, 
they can be divided in two different classes, on the basis of the behaviour of the frequency at early 
times and the value of the Bogoljubov parameters. 

The first class is represented by those dispersion relations for which there is no violation of the 
WKB condition at early times, but only at very late times when the dispersion has the linear form 
uj ~ k as in the standard case. A good example within this class is the Corley-Jacobson relation with 
positive coefficients [22] . which was studied in [21] , In this case, the WKB approximation (well-suited 
for these relations) shows that the corrections to the spectrum and bispectrum are very suppressed. 

The second class of dispersion relations consists instead of those violating the WKB (adiabatic) 
condition at early times. This scenario is more promising: first, interference terms will generally appear 
in the bispectrum due to the presence of both positive and negative frequency solutions in consequence 
of the early WKB violation. In this case, one could therefore expect large enhancement factors for 
some configurations because of interference and phase cancellation. Second, the Bogoljubov coefficients 
indicating particle creation are actually quite independent of any ratio of scales of the form as we 
will show, and are proportional to the parameter signalling the WKB violation and the interval of 
time during which the WKB approximation is not good. The only constraint is due to the request of 
small backreaction [llj . 

In the following, we will study the corrections to the bispectrum for dispersion relations with the 
behaviour represented in figure UJ where the WKB condition is violated at early times. The times 
"Hit Vn-i signalling the transition between violated /satisfied WKB condition are individuated by the 
violation of (some of) the conditions in equations (|14p . (|15p . and are close to the turning points 

5 The slow-roll parameter s = —j^i^f should not be confused with e = 

6 It might happen that with an appropriate change of coordinates the WKB approximation could be improved also 
in the region where it would not be good when using the conformal time variable r\ or the rescaled variable y = e kr/, see 
for example what shown in [2H for the case of standard dispersion relations. However, we will consider the general case 
and keep using rj, y and speaking of WKB violation when the WKB conditions are not met in these coordinates. 
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Vii V2, V3 m figure [TJ They are inversely proportional to the wavenumber k, being given by a condition 

(k) 

on y = ekrj, and hence, more correctly, we will indicate them with the notation ry- when necessary to 
avoid confusion. It is important that at very late times the WKB condition (adiabaticity) is restored, 
because in that case we can unambiguously determine a vacuum state and fix initial data. We will also 
briefly comment on the case of dispersion relations with no WKB violation at early times, generalizing 
the results found in |21j . 

Dispersion relations that are more complicated than the one in figure [1] can be imagined, with 
several period of WKB violation and a more varied behaviour, but the analysis in those cases amounts 
to repeating the one we do for the case in figure [1] in the various regions of WKB validity or violation. 
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Figure 1: Dispersion relation with violation of WKB at early times. 



3.1 Solving the field equations 

In the case of a dispersion relation like the one in figure [IJ the time axis can be divided in four regions, 
and the solution written as follows 

, B 1 U 1 ( V ,k)+B 2 U 2 (v,k) 
jt{ri) = < 

a^ui(rj, k) + P%u 2 (ri, k) 
DxVi(v,k) + D 2 V 2 {v,k) 

The coefficients a?, /3^, -61,2, -^1,2 are determined by imposing continuity for the function and its first 
derivative, is determined by the choice of boundary conditions for r\ — > —00. 

We also must impose the Wronskian condition W{f, /*} = —i to have the standard commutation 
relations in the quantum theory. If we impose it at a certain time, the continuity of the function and 
its first derivative ensures that it will be always satisfied because of Abel's theorem and equation Q. 

The partial solutions listed in equation (|24p will be presented in the next sections; in particular, 
we will use the WKB method to obtain u\ % and choose the adiabatic vacuum for the field. 

3.1.1 The WKB solution of the field equations in regions I and III 

As we have shown, the WKB method has lead us to the equation (|13p with q(y) given by (|2ip . which 
can be solved in sequence comparing powers of e. In the regions I, III the correct distinguished limit 



I : 7] < rj{ 

II : 
III 



CO 



0) / ^ (k) 
T)i < r] < 7]x 



(*0 ^ ^ (k) 



(24) 



IV : rS ] < rj 
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is 5 ~ e, leading to the sequence 



s 2 + F{y) 2 = 
2s si + s = 
sf + 2s S2 + s'i 
2sis 2 + 2s S3 + S2 



2 

72 



(25) 



To take into account the — term in the field equation, we need to include up to the order 5 2 in 
the series expansion for the WKB solution. The Wronskian condition then forces us to consider also 
the order 5 3 . In this way we obtain the positive- and negative-frequency WKB solutions^ u\ 2 



ui(Vk) 



■£n(yjb) 



y/2kU{y k ) ' 



uziVk) = u*{y k ) 



■|n(v*) 



y/2kU(y k ) 



U(y' k )dy' k , U(y k ) = S + e 2 S 2 = F + e 5 



Ip2 



+ 3 



8F 3 



(26) 



(27) 



With the above notation for £l(y k ) we intend the primitive of U evaluated at y k = ekr], since we 
absorb the phases resulting from the lower limit of integration in the coefficients o^, ap, /3r. The 
notation y k keeps track of the /c-dependence of the definition of y and will be useful in the following. 
We have also normalized the solution with a factor yk because in this way the Wronskian constraint 



yV{u\,U2~} = —i in the //-coordinate is solved if ?g is a pure phase, and 

The WKB approximation is good when all the conditions (|14p . fjl5j> are satisfied. In particular, 
since we truncated the series at the order <5 3 , we must require, among the other conditions, 



< 1. 



(28) 



Observe that, by direct derivation, our WKB solutions are found to precisely satisfy the equation 
(using the coordinate rj) 

< 2 + [V(V, k) 2 - Q( V , k)) ui, 2 = (29) 

where 



// '2 



1 V co" 



J 2 U" 



+ -3-r-ir - — +3 



U 



'2 



-e A 2s {y k )sA{yk) , (30) 



4w 2 8w 3 wy 2 J 2to 4uj 2 2U 4U 2 v^vt 
therefore the condition (|28[> corresponds to | | <C 1. The times rj u ,rj m are then such that | ^% | ^ 1. 



3.1.2 The solution of the field equations in regions II and IV 

The WKB approximation is not good in the regions II and IV |§ In previous works, then, it has then 

often been taken the drastic approximation that V(r], k) 2 ~ in those regions, so that the solution of 
equation ([9]) reduced to a growing and decaying mode depending on the scale factor. For the purpose 
of studying the bispectrum, we employ more refined and accurate approximations. 

7 From the ansatz JUJ of the WKB solution we actually obtain e ± '" (vk) ^^ log(So) " V If i nst ead than 
e ±|n(«)-3 iog(s +e 2 s 2 ) j but the resumme d f orm i s better for the Wronskian condition. Here, recall U = So + e 2 S 2 - 
8 See however footnote [6] 
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Indeed, in region IV, when rj > r/ m , the dispersion relation is in the linear regime 

w 2 (r? > r] m , k) = k 2 , 

r(l,2). 



and the solution is given in terms of the Hankel functions V± j2 = \J—f}Bv ' (—kij), v ~ |. 
for the continuity of the function and its first derivative at r] < r/ m , we find that 



2 k 



(31) 
By asking 



(32) 



When instead % < rj < ?y n , in region II, we indicate with Ui t2 {k,rj) the solutions obtained not 
neglecting the uJ 2 term in the potential V 2 in equation Q. Actually, we will find that some of the 
important features of the bispectrum can be found without knowing the details of the field solution, 
but just using the fact that small backreaction requires having a short interval [rfi, as we will show 
in section [5j In this sense, those results will be very robust, being independent of the approximations 
chosen when solving the field equations. 

When needed, however, a rather accurate approximation in region II is obtained observing that, 
given figure [H the physical frequency has a local minimum ujq at r] = rj m , and therefore 



2 2 

a ui , ~ 

phys 



V 2 H 2 



ul +w w phya (r/ m )(7/ -?7 m )' 



(33) 



This approximation is well-justified also because, as we said, backreaction constraints the interval 
[?7i! fjji] to be very small, and therefore r\ ~ r] m for rj, r\ m G [r^, see section[5]for a proof. Substituting 
(|33p in Q, the solution is readily found: 



fArj) = Bi W(iKrj m ,a,2iKr]) + B 2 W(-iKr) m ,cr,-2iKri) 



K " H 

where W(a, b, z) is the Whittacker function. 



a = 



^9H 2 - 4lo 2 - AH 2 K 2 rj^ 
2H 



(34) 



3.1.3 Initial conditions, the Bogoljubov coefficients a^, (3^, Bi, B 2 and particle creation. 

The result regarding the coefficients a^,(5^,B\,B 2 in (|24p is robust, because it is independent of 
whatever assumption is made about the behaviour of cj, in particular, of the form of u assumed in 
equation (|33|) . The only information we need is that the interval [/ft, %] is very small because of 
backreaction (see section [5]) so that 

rji ~ % 



A 



rji 



< 1 



and that the WKB solutions satisfy the equation (|29|) . while the partial solutions U12 satisfy 

<2 + V(r), k) 2 U h2 = rjeh,Vn}- 
By asking for the continuity of the solution and its first derivative, we obtain in full generality 



B> 



a j 



W{Ut,U 2 } 
W{B 1 U 1 + B 2 U 2 ,u 2 } 



m 



W{U!,U 2 } 



Bo 



W{Ui,U 2 } 
W{B x U x ^B 2 U 2 ,u x } 



m 



W{ux,u 2 } 



(35) 
(36) 

(37) 
(38) 
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where W is the Wronskian. By expanding for small A around we obtain (again, in full generality) 



V 2 



A+0(A 2 



in 



St' 



V 2 



A+0(A 2 ) e (39) 



but to the 



We see that at leading order /3g is not proportional to small ratios of scales such as 
parameter ^ | signalling the WKB violation (see equation (|30p ) and the ratio A. The only constraint 
to the value of &c is therefore due to the request of small backreaction. 

The Bogoljubov coefhcients also depend on the choice of initial conditions, that is on the pure 
phase The bispectrum will be affected by this choice as well, differently from the spectrum, where 
only the modulus square of the mode functions enter the computation and therefore |cr| = 1. In the 
following, we will choose Cf = 1. 



4 Bispectrum and modifications to the Bunch-Davies result 

The three-point function for the comoving curvature perturbation can be computed in the ip ^ 
gauge and then transformed in the £ 7^ gauge, using equation ([3|) at leading order. In the following, 
when we write a formula in terms of £ 7^ it must therefore be intended that the computation is 
performed in the ip 7^ gauge and then transformed in the £ 7^ one, where the quantities relevant 
for observation are more easily read. 

In the in-in formalism, one obtains at leading order [2] 

(C(xi)C(s2)C(*3)) = -2Re ( f d7 ? / i^in|C(xi)C(x 2 )C(x 3 )%)(V)l^m)') (40) 



> hi 



where Hrj\ is the interaction Hamiltonian, while rji n , \ipi n ) are the initial conformal time and state 
(vacuum). The standard result in slow-roll inflation is obtained by choosing \ipi a ) to be the Bunch- 
Davies vacuum and r/; n = — oo [2]. 

The basic correlator in perturbation theory is given by the Whightman function 

We write the bispectrum as |25| 

(Cfc (vKk 2 (V)(k 3 (V)) = F(h,k 2 , k 3 ,v) = (2vr) 3 5(^ h)B$iM, 4 v) ■ (42) 

i 

Scale-invar iance requires the function B to be homogeneous of degree —6 and rotational invariance 
makes it a function of only two variables, which can be taken as x 2 = |^ and X3 = The conservation 
of momentum, with the triangle inequality, forces x 2 + £3 > 1. Because of the symmetry in x 2 and 
X3, it can be further assumed that X3 < x 2 . 

Typical configurations in the standard slow-roll/chaotic models are the local one, where x% ~ 
and x 2 — 1, and the equilateral, where x 2 — X3 ~ 1. The standard inflationary model with Bunch- 
Davies vacuum state tends to produce non-Gaussianities when the modes cross the horizon, leading 
to a three-point function maximized on the equilateral configuration. If instead the non-Gaussianites 
receive a contribution also during the superhorizon evolution, as for example in the curvaton scenario, 
they will tend to be maximum on a more local type of configurations |25j . 

In the next two sections we compute the three-point function for two cases of interactions: a cubic 
coupling without higher power of derivatives acting on the fields, and a cubic interaction with higher 
derivatives. Some of the details of the computation apply to both cases, and we discuss them at length 
within the former. We then use that knowledge to analyze the latter one, which will turn out to be 
the most interesting for observations. 
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4.1 Cubic scalar interactions 



The interacting Hamiltonian at cubic order in the perturbations can be obtained by expanding the 
Einstein-Hilbert action on the quasi de Sitter background [2j. The contributions to the correlator 
(|40p include a connected and some disconnected parts depending on how we define the curvature 
perturbation (that is, if we use non-linear field redefinitions). We will be using a redefinition leading 
to the simplest form of the cubic interaction [2J: 

%) = -/^ 3 - 3 (#) 4 ^C' 2 5-V, (43) 

Planck 

and will be interested in the connected three-point function. In the following, we indicate with the 
terms "three-point function" and "bispectrum" their connected parta^l. 

Using equation (03]), the three-point function (|40p becomes, in momentum space, 

faKfc (r))(n 3 (V)) = 2Re (-*(2vr) 3 ^ (jf) ^ f d ^^ff~ II d v' G k t (V, r/)+permutations 

\ i \ / Planck^^in 3 j — 1 / 



where r/ is a late time when all modes ki are outside the horizon after the modified dispersion relation 
has become effectively the standard linear one. 

Evaluated using equations (|24p . (j32|) at the late time \rj\ <C 1, when all three mode functions have 
exited from the horizon in region IV, the Green function reads 

To compute the bispectrum, one needs to divide up the interval of time integration in equation 
(|44p into the different regions of validity of the piecewise solution (|24[) and compute the various terms. 
The number of terms is elevated, but many of them are identical up to permutation of the external 
momenta and therefore we can reduce them to some common classes which we now study. By the 
shortcut expression "being in region X" referred to a Green function depending on fcj, we will mean 
that for those values of rf,ki the function /^(V) entering equation (|45p is the partial solution of the 
field equation valid in region X, as shown in equation (|24p . 

4.1.1 All Green functions in region III 

We start from the analysis of the contribution to the bispectrum ()44p when all the three Green 
functions are in region III, which, as we will show, leads to the largest enhancements. In that 
region, the function in (|45p has the form 

drfGfarf) = fyGjfW) + d V 'Gf\v,v') , (46) 

where we have distinguished the positive- and negative- frequency branches ui^(k,r]) of the solution 
in region III, see equation (|24p . 

There are therefore two possible contributions to the three-point function: one where the three 
Whightman function are all in the same positive or negative branch of the solutions, the other where 
one of them is in the opposite branch. It is convenient to concentrate first on the latter case. 

9 Therefore, our results can be confronted with the observations only after adding the disconnected contribution, which 
in any case do not have enhancements and are subleading in slow-roll parameters, see [2]. 
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The leading contribution in powers of f3^ to the standard slow-roll result is then 

H 2 H 



3,y3 



-2Rc 



1(2^5(^2^)- 



i 2 M 2 

,t JW Planck J Vk 



f im 1 e^H 3 T-r 



fc ?nti(2^ 



Mi 



x e 



+ permutations /i, j G {1, 2, 3} . 



(47) 



The Green function depending on kj is in its negative-frequency branch, and those depending on 
kh, h 7^ j are in their positive-frequency one. The limits of integration r)' u ,r)' m are given, respectively, 



by the largest and the smallest among rfr^ and 77 



„(*0 



1,2,3. We have also defined 



n ^ 1 ( 1 U' .__ W\ 



(48) 



and the functions U(rf,k) = U{yk{rf)) and 0(r/, k) = £l(yi~(r)')) are given in equation (|2"T|). "H is the 
Hubble rate in conformal time, M„, , is the reduced Planck mass. 

' Planck 

We rewrite (j4"T|) as 

F 3 ^(fci, k 2 , h) = A SF 3 ^(k 1 ,k 2 , h) (49) 
where we factorize the amplitude, including the overall scale dependence, 



A = 4(2vr) 3 5(^fc 



h 2 h 2 h 2 



Planck IlLl(2A?)Y** 



1 3 



i=l 



while 



8F3,p(h, k 2,k3) = ^Re 



0i 



drthX[l*$ h ,Tt)7(kj,i{)e i 



(50) 



(51) 



The amplitude A is precisely the standard slow-roll result [2], hence a is the putative enhance- 
ment factor, to which we now turn. Recall that the modified dispersion relation implies that the 
frequency oo(j], k) depends on k, t] and the scale of new physics A as in equation (|19j) . We then change 
the variable rf in the integral in equation (|5ip to 



- H i- < 



max(&i, k 2 , ^3 



so that equation ([51]) reads 

8F 3t p{xi = l,x 2 ,x 3 ) = ^Re 



P% A 

-YH Xt 



urn 



(52) 



(53) 



where we have assumed, for instance, k\ = /c max and recovered the usual bispectrum variables x% = , ' 
introduced at the end of section [2J To write the formulas in a more symmetric way, here and in the 
following, we also indicate formally the ratio x\ = , kl , although it is equal to 1 because of our 

""max 

assumption. 

We have also defined 



Sp(x!,x 2 ,X3,y) = ^2 x h ) - tt(y, xj) = S /3({x},y) + -^S 2 p({x}, y) 



S (2)p({x},y) =^2,S 0(2 )(x h ,y) - 5 , (2)(^,y). 



(54) 



(55) 
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where the dependence of Q on Xi,y is due to the fact that yki(v') = ki er }' = x iV? see equation ([2] 
and So, S2 are the quantities appearing in the second equation in (|27|). 

The limits of integration in (|53p are computed as the values for which the WKB conditions are 



violated. In general, one finds that y n = x^max^n m ~ ~~ 1> since the corrections to the linear 
dispersion relation at the time r/ n must be quite important in order to drive the frequency close to 
the turning point (see figure [TJ, and therefore looking at f)19|) it must be S2^l ^ 1. The upper limit is 

instead y m = f k^J^ ~ -f 



mm 



Evidently, the contribution (|53p is sizable if the interval of integration is sufficiently large. We are 
therefore interested in this case, hence we will consider x. > lOe, which however for e reasonably 

fc min 

small leaves plenty of room for all the values of the wavenumber relevant for the CMBR observations 
(consider for example that already for the supersymmetric GUT scale, it is e < 10~ 3 if H ~ 10 -5 Mp). 

The integral in (|53p is a typical Fourier integral which can be well approximated by the technique 
of stationary phase, since S> 1 [26] . As it is known, its approximated solution depends on the critical 
points of Sop({x},y) as a function of y. Their presence and nature depend on the configuration of the 
external momenta k\ 2,3, which act as parameters in the function. 

We call stationary point of order n — 1 an interior point such that 

d n y S w ({x},y)\ y=y , = S^({x},y,) 0, d™S p({x},y)\ y=y * =0 Vm<n. (56) 
The leading order solution to the integral is then [26] 
- if n e 2N 

= £ (iV~" ,iff r ; ite[nr(x fc ,i/07( g3 ^o^^ < ^ (W ' y -^^^ )({ ^ ) ] 



('^ ^, ')" 



(57) 



- if n G 2N + 1 > 1 



- if n = 10 

^ = ^E g( i)"\ i m [n^fa.^^)^/ 4sMW - y) ]i;: (59) 

i ^0/3 (l x )'y) h£j 



Observe that j(x, y) does not go to zero if the WKB approximation is valid. 

We recognize the appearance of enhancement factors proportional to (-jj) 1 "in the presence of 
stationary points. We postpone the detailed analysis and comments regarding this to section 14.1.61 

10 In fact, if \S' p({x}, j/ni/n)| < \f^{ 013 ~ ^' ym/E ^ ) 3 ) the correct approximation must take into account the higher 
order terms in the expansion of Sop. A better result including the second order is then given by (|57[) for n = 2, evaluated 
at Vm/Ji, divided by -2, and further multiplied by e~ l ~ slenS o.ft tW-^ni' m tjjg case Q f y m jf necessary, one can also go to 
higher orders. In these cases the boundary point can be called a nearly critical point. 
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4.1.2 All Green functions in region III in their positive-energy branches and the case of 
dispersion relations not violating the WKB condition at early times 

The contribution to the bispectrum when all the three Green functions in formula (|44p are in their 
positive-energy branch can be obtained by the replacements 

-+ a| ~ 1 , -iSltf, kj) -> iSl(r/, kj) (60) 
in equations (|47p and successive, which in particular amount to the replacement 

Sfi({x},y) ->■ 5 a ({a;}, y) = 5 0a ({x}, y) + 5 2q (M, 2/) = ^ S (xj,y) + ^ S 2 (xj,y) (61) 

3 3 

in equation (|53p . If we now investigate the presence of stationary points for the function So a ({x},y) 
as we did before for Sop, we find that there are none. In fact, let us check the vanishing of the first 
derivative, which is a necessary condition to have a stationary point of any order: 

Soa({x},y*) = w(xi,y*) + u(x 2 ,y*) + u(x 3 ,y*) = 0. (62) 

This condition is never satisfied for non-trivial configurations, as the quantities u(xi,y*) are always 
positive. The approximated solution for the contribution to the three-point function is then 

6h« = \T. g(1) *\ I m [[I^fa-^^)^ Sa(W ' y) ]e (63) 

The case of a dispersion relations not violating the WKB condition at early times follows this same 
pattern with /3r = to start with, while the time integral in the analogous of equation (|53p is extended 
to minus infinity^. The necessary condition for the presence of stationary points is therefore still given 
by (|62p . which shows that none appears and therefore there is no enhancement, in agreement with 
what found in [21] . 

4.1.3 Green functions in region IV 

By computing the relevant Green functions, it is seen that the contribution to the bispectrum when 
(some of) the Green functions are in region IV can be obtained by replacing the functions y(k, rj) and 
U(k, rj) in equation (|47[) respectively with i and 1 for the Green function(s) in region IV. However, 
since we imagine computing the bispectrum shortly after the exit from horizon of the modes, because 
in that case the error involved in using the relation Q at leading order is minimized, this contribution 
will not be very large. 

4.1.4 One (or more) Green function(s) in region II, the others in III 

We now briefly study the contribution to the bispectrum when one (or more) of the three Green 
functions in its formula are in region II. For definiteness, we consider the case where one Green 
function, say the one depending on k$, is in region II and the others are in region III. It follows that 

ki, 2 < h . (64) 

The other configurations are readily obtained by permutation. 

11 The time integration path must be chosen such that the oscillating piece in the exponent of the integrand becomes 
exponentially decreasing. This corresponds to taking the vacuum of the interacting theory. 
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The generic contribution to the standard bispectrum can then be written as 

(fc 3 ) 

x 7 W (A?! , r ? ')7 ( * ) G k:i (r), rf) e^^' M)±&n(v' te) ) + permu tations. 



where the signs in the phase of the integrand are plus/minus if the contribution involves (a//3)* , 

fcl,2 

and similarly we have 7*/7 in presence of (a//3)*. . 

The three-point function includes now a Green function of the form, see equation (|45p . 

„ „ , lfie-iH j B\Ul(,j ,k) + BjW 2 (rj ,k) \ 

W** > = ^3 ( W) ) ■ (66) 

By using the results in (|37p . we obtain 



n 2 fa) ^fa^j^O + KMVi) U^Ml^ri') 
+ a(r/)W{UUVx),UZ(rh)} 

where [ ] indicates antisymmetrization. 

The time integral in equation (|44|) is over the interval [Vi k3 \ Vn^ ] = [Hi^ ; ^i*" 3 ' (1 — A <fe3) )] , thus 
we change variables as 

rj' = 7,^(1 -A (ka) v) (68) 
and obtain, after expanding for A <C 1 and recalling that V(r]i, k^) 2 ~ Ti 2 , 

dG k3 (v ~ 0,1,) ^ e^^4^^1(l + *A (fc3) „) . (69) 

We now write 

^3,B = ^4 5^3^ + permutations , (70) 

where A is defined in (|50|) . We also change the variable y' k . in (|27p to £ = A^.fa)'' so that the 
equation is more naturally written in terms of the variable v defined in (j68[) as 

— U(y ki ) -> -? ?I A / kiU(Z,ki)dt = -r}iA hU(xi, h) dt = - Vl A n(«, A*) , (71) 
and equations (f65|) . (f70|) give 



5F3.B = Re 



A (fe3) ( fc 3) (*3) 



2 6 



&3 ^3 



(72) 



where 

K B fe,fc2^) = T%fel)T%fe), =F**(a/0)c • (73) 
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We see that now the asymptotic stationary phase approximation of the Fourier integral is not accurate, 
as there is no large factor in the phase of the integrand (since A (fe3> <C 1). Therefore, also when the 
combination of signs leads to interference and phase cancellation, there will appear no enhancement 
factor. Observe also that the correlator is proportional to the small factor A (,=3) <C 1. 

It follows that, overall, this contribution is suppressed compared to the one we discussed in section 
14.1.11 It is straightforward to realize that when more of the Green functions entering equation ([44 p 
are in region II, the contribution is also suppressed. This result is very robust, because it is not based 
on the specific form of the solutions k): we have in fact only used that A (,=3> <C 1 and the general 
form ()37p of the coefficients B\ 2 which is valid for all solutions. 

4.1.5 One (or more) Green functions in region I 

The last kind of contributions that we are left to analyze is the one where one or more of the Green 
functions in the formula for the bispectrum are in region I. There can be various possibilities, depending 
on what region the other Green functions belong to at those times, given V(r],k). If (some of) the 
other Green functions are in region II, we expect the contribution to be suppressed in a way similar 
to that discussed in section 14.1.41 Recall that in region I, the WKB solution has only the positive- 
energy branch u\ because we chose the adiabatic vacuum, see (|24p. therefore no interference is possible 
between functions in region I only. If instead at least one of the Green functions is in region III, we 
could conceive the presence of interference terms possibly leading to enhancement factors 12 !. 

Let us consider this case in more details, imagining for definiteness that the Green function in 
region I is the one depending on ks, while the others are in region III. The form of the contribution 
to the bispectrum is the same as that in equation ()47p . with the difference that now w(ry, ^3), 0(77, ^3) 
and 7(^3,77) feel, through equation (flllj) . effects that are not suppressed in region I. 

Following the derivation in section 14.1.11 in order to have the enhancements there must exist 
momenta configurations allowing the existence of stationary points for the function given by the 
sum of phases of the WKB solutions multiplied by the large factor jj . In the situation we are discussing 
here, it however appears more difficult to meet this requirement. Indeed, let us consider just the first 
and necessary condition for having stationary points of order at least 1: the equation is, say, 



We expect that satisfying this equation for non-trivial configurations would be more difficult than in 
the case of section 14.1.11 because the frequencies are now in different parts of their curves: u{x\^) in 
region III, while uj(x%) in region I. Therefore, we expect not to find enhancements. 

However, since this obstacle to having large enhancements strictly depends on the form of the 
frequencies uj, there could be cases where the conditions for the stationary points can be satisfied. As 
this is strongly model-dependent, we will not deal with this point further, and leave it to be addressed 
case by case in the models of interest. 

4.1.6 Comments on the features of the corrections to the slow-roll bispectrum from 
Lorentz and WKB violation at early times 

It is convenient to pause for a moment and comment in general terms the results found so far. The 
most notable features of the bispectrum from cubic interactions in the case of Lorentz and WKB 
violation at early times emerge from the equations (|57|) . (|58p . (|59p and (|63p . as we have shown in 
sections 14.1.31 14.1.41 14.1.51 that the other contributions are suppressed. 

12 For this, at least one of the Green functions in region III must be in the negative frequency branch of the solution, 
otherwise there will be no enhancement, as discussed in section ELl.il 




(74) 
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The contribution given by equation (|63p . multiplied by the factor in (|50p . matches the standard 
result, from which it deviates only for the presence of tiny superimposed oscillations. The same 
qualitative behaviour of the bispectrum occurs for dispersion relations not violating the WKB condition 
at early times, which therefore do not show strong deviations from the standard slow-roll results. 

The corrections to the standard result for dispersion relations that do violate the WKB condition 
at early times are instead quite different. They present a modulation in wavenumber space, descending 
from the oscillatory behaviour of ([57]) . ([58]) . ([59]) . but in addition the contribution of certain configu- 
rations can also be enhanced, depending on the form of the dispersion relation. These configurations 
are individuated by the system (|5(jj) . and can therefore be found by the knowledge of the dispersion 
relational. Compared to the Bunch-Davies slow-roll result, they lead to the enhancement factors 

ga = ^~WAY~' 1 |^|, ,75) 



BD 



which can be large for n > 1. 

Let us compare this result with the one found in [3] and anticipated in [9], based on the modified 
vacuum approach, which models the transplanckian physics imposing a cutoff on the theory in mo- 
mentum space. The result in that case also showed oscillations in the bispectrum and enhancements 
of the order of -jj, but only for the enfolded configuration. We see that in general the enhancements 
in the case of modified dispersion relation have a smaller magnitude: it would be the same if there 
were configurations {x} e nh> such that there existed a stationary point y* with the property that all 
derivatives of Sop({x} eTi h, y) at were zero. The enhancement factor would then be given by the 
limiting behaviour of (|75p for n — > oo. In fact, this occurs for example if 

Sop{{x},y) = Sht{{x},v) = (^2%h - Xj){y - Vc), Vc = —m&x(ki,k 2 ,k 3 )r] c , (76) 

which is indeed the phase of the integrand in the bispectrum formula obtained via the modified vacuum 
approach, where rj c is the time cutoff corresponding to the appearance of new physics [3]. In this case, 
all derivatives of Sfrriik}^ x) are zero except the first, which is zero for the enfolded configurations 
given by Sfe^j ~ kj = 1 1 I. However, from the point of view of modified dispersion relations, where 
uj 7^ k, it appears that such an eventuality would be rather peculiar. 

The final magnitude of the enhancement depends on the value of the Bogoljubov parameter 
In section [5J we will show how the constraints from backreaction affect it and reduce the magnitude 
of the overall enhancement factor ([75]) . However, by looking at equation ([53]) . we spot that larger 



enhancements would appear for interactions that scale with higher powers of We turn now to 
one such example. 

4.2 Higher derivative interactions 

An example of interaction that scales with larger powers of than the cubic coupling fi3"[) is [31127] 



£ (J) = v/Zdet^^V^) 2 ) 2 . (77) 



Once expanded in perturbations, this term leads to one contribution quadratic in the fields and 
to another which is cubic. The first one has the effect of modifying the sound speed of the modes 



13 Also near-to-critical points lead to enhancements, which are however smaller than those of critical points, see footnote 
EH 

14 In this case we cannot talk of stationary points, since they are not isolated: the derivative is zero over the whole 
range of y. 
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in the field equation ([9]). We neglect this and focus instead on the three-point function generated by 
the cubic contribution to the Lagrangian and, consequently, to the interacting Hamiltonian. It can be 
found that the latter is [3] 



%) = -j is? - (^) 2 ) =-j d^a^C'(C' 2 - (S.C) 2 ), (78) 

From equation (|4(jp . we obtain at leading order the correlator 



F(k!,k 2 ,k 3 ) = 2Re 



(2n) 3 5^ (J2 h) ^4 J d V 'a(d v G kl (r,, 77')^ (77, v')d v G k3 ( V , v')+ 



+ h ■ A; 2 G fel (r?, rf)Gk 2 (rj, rj f )d v G k3 (77,77')) + permutations = A 5F d , (79) 

where A is defined in equation (|50p . 

From the analysis in section |4~T| we have learnt that the largest correction to the bispectrum occurs 
when all the Green functions are in region III and one among them is in the opposite frequency branch 
with respect to the others. We therefore focus on that case. At leading order in /3^, 

\ hft X j + 1 X j+ 2^U{x j + 1 )U{x j+2 ) X j+1 Xj y/U{X j+1 )U[Xj) 



+2 



x j+2 ■ Xjj*(x j+1 ,y) 



Xj+2 Xj yJU{xj + 2)U{xj) 



(80) 



where we have changed variables as in equation ([52]) , written X{ = , i and used 7,673 defined in 

"'max 

equations (1151) . (1511) . while C = ^ a^l . We have not indicated the y-dependence of 

U(xi,y) to avoid cluttering the formula. In the expression for P({xh}, xj, y), indices are defined 
modulo 3. 

Now, we need to study the presence of stationary points for the function Sop, as done in section 
14.1.11 We need also to check the behaviour of P({xh},Xj, y). In fact, it could go to zero or be very 
suppressed for the configurations for which there exist stationary points of Sop. This is what happens 
in the modified vacuum case, where the highest power of - in this bispectrum contribution is exactly 
zero on the enhanced folded configurations, and therefore the final enhancements is reduced [3]- 

It is straightforward to observe that in the case of modified dispersion relations, P({xh},Xj,y) does 
not generally go to zero nor it is strongly suppressed. Let us consider the case of a stationary/boundary 
point at the value y = y* ^ 0. Since |y*| < |y E | ~ 1 we can expand^! P in powers of y*: 

l(xi,y*) ~ i + iciXiy* + — ) + . . . (81) 

z Xiy* 

where Cj = d y U(xi,y)\ y= y t ~ 0(1). It follows thatr^l 

iP({x h },Xj,y) ~ k + + k oa y* + k 10 e + ac u ey* + . . . (82) 

y* 



15 The expansion is not well justified if y* = j/n, but even in that case one can check that there is no cancellation of 
P({x h }, Xj,y*) in general. In the following, recall also that < 1 for all y» £ [ya, ym\. 

16 In equation ()82p . indices with the letter h are defined modulo 2, indices with the letters j or i or m are defined 
modulo 3. 
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XlX2 X 3 - - - - ■ <-*i=i »+ 2 

3 

K a,i = S^^c^i — (^^Cfta;h^{;j +1 a; /i+2 j — CjXjX^ +1 x. +2 -^) + (^^ChX h x h+l x h+2 — CjXjXj +1 Xj +2 ) 

i=1 hjtj h^j 

+J2 x h{-c h +i + Cj) + x^c h + ^ 9 h h (x 2 h (x h+1 - x 3 ) - x 3 h+1 + xf) + 3 ( x "jJ2 Xh ~ H x h) 

h*j h#j h& XlX2X3 4C1X2X3 ^ ^ 

1 3 

K i,o = 3i(y Cj — Cj) — ]chXh(x h+1 — Xj ) + CjXj(Xj +1 + Xj +2 )) — ~^^ J { c i x {h+l x h+2} CiX i X h+l X h+2) 

h=^j h=£j i=l 

K i,i = Zi-Cy^ChXhjxh+i - Xj) + CjXj(xj + i + x j+2 )) + Giicj+iCj^Xj+iXj^xJ 1 - CjCj + iXjXj + ixJ+ 2 

3 3 
/V^ , CiCmX {i X m} v 

- CjC ]+2 X 3 X ]+2 X jJrl ) + -(2_^XiC l+1 Ci + 2 + 2^ ) . 

^ -, . X1X2X3 

where { } indicates symmetrization. 

The coefficient k matches the result of [3] as expected, and it is indeed zero for the enfolded con- 
figuration. In general, however, in the case of modified dispersion relation the enhanced configuration 
need not be the folded one, and therefore k 7^ 0. Even in the case of the enfolded configuration, the 
other coefficients of the expansion will not be zero in general and need not be utterly small. 

Therefore, for a stationary point of order n — 1, we obtain [26] 

- if n E 2N 

6h, p = -£c (A) 3 " " J£|L _ Re [P{{XhH} , Xj , Vm) p % jtom^^&i*}^ 

(83) 

- if n £ 2N + 1 > 1 

= -TP (£)'"* J^L^Mm^MH^"^ ■ ^ 

\ n\ ) 

We conclude that, in the case of modified dispersion relation, the enhancement for higher derivative 
interactions, when present, is actually larger than the one obtained via the modified vacuum approach 
to transplanckian physics, if 

A \ 1-1 

jj) y* 2 >i, (85) 

as can be seen comparing with equations (3.31), (3.32) in [3]. We comment on the likelihood of this 
at the end of the paper. Finally, 

- if there is no stationary point, n = 1 m 



17 If \Sop({x}, 2/m/n)| < ^/x( " ° g ~ ^' Vm/n ^ ) 2 1 the correct approximation must take into account the higher order 
terms in the expansion of Sop- In this case, a better result including the second order is then given by ()83|) for the case 
n = 2, evaluated at ym/n, an d then divided by -2 and further multiplied by e _l ~ slgnS o,/3 (M-fin) m ^he cage D f y m jf 
necessary, one can also go to higher orders. However, observe that ym is such that certainly (-^) 2 yni ~ lj an d therefore 
there is truly no enhancement at that boundary point. 
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Schematically, the enhancement is given by 



Fed 



3 



s$({k},y*)\ 



-\Pz 



(87) 



5 Backreaction 

The energy density produced in connection with particle creation backreacts on the background. When 
interactions are turned on, this energy can be divided in a "free" and an "interaction" parts. 

The study of the backreaction in the case of modified dispersion relations has so far only dealt with 
the free part, as in |12]I28] . The interaction one has instead been studied mostly within the approach 
of modified initial vacuum in [29] and especially in [3]. We will show that the general results found 
in [3] are similar to those we find for the case of modified dispersion relations. 

We start with a review of the features of the free part of the produced energy density. We consider 
the energy generated by the particle creation in region III, which is the most relevant for our results. 
The formula for the average energy density can be obtained computing the energy-momentum tensor 
of the inflaton from the typical Lorentz-breaking action we report in appendix IA.1I [12 : 
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Using the third line of the solution (|24p . valid in region III, one obtains 
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(89) 



where the momentum integral extends over values for which the solution of region III is valid. Here, 
g = \^+ikU + U ~ iuj(k,rj), being in region III, so that for modified dispersion relations the terms 
with exponentials cancel at leading order even without doing the integration (averaging) over the 
oscillating exponential, which, as argued in [3j, would damp those contributions anyway. 
The formula for the "free" energy density is therefore reduced to 



(0|p|0) 
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W phvs(P). 



(90) 



in terms of the physical momentum V = \ an d frequency u phys (p) = = pF(—^) , see equation (|19|) . 
We now substitute the expression ([39]) for the Bogoljubov parameter and change variables as 
= y{rj), so that the integral ([90]) . after discarding the zero-point contribution as usuaf^l. yields 



2 A 4 



(0|p|0) ~ A^A' 



(91) 



This has to be compared to the background energy density pb = 3H 2 M 2 , yielding the constraint 



H M P1 . 

A _s .Planck 

A 2 



(92) 



18 In a gravity theory this is a non-trivial step, but we uniform ourselves to what usually done in the literature on the 
subject [121128] . 
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to avoid issues with the backreaction. It ought to be required that the minimum ojq of u h s at early 
times (see figure [1]) be different from zero, otherwise, as discussed in [28^30] . there would be WKB 
violation even at present times and the constraint would involve H today i becoming very constrictive. 
If instead ujq > 0, it can be shown that there is no further backreaction issue past a certain time after 
the end of inflation, because the Hubble rate decreases and the WKB violation does not occur any 
more after the time when H WQ ~ ujq. In our estimates in this section, we make the assumption that 
Huj ~ H at inflation for simplicity, but the general case can certainly be considered. 

We also have to preserve the slow-roll conditions. We define the standard slow-roll parameters and 
link them to the energy and pressure density via the equations 



dH 

lit 



-eH 1 



1 



2M 2 

Planck 



(p + p) 



d 2 H 
dt 2 



1 



2M 2 

Planck 



ft " mp + P) 



where we have called fi = r] s — e, and rj s is the slow-roll //-parameter. 
Using the formula [12] 
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(94) 



we obtain (p) ~ A A and ( 



>d£\ 
dt I 



A 2 HA 4 , from which the bounds 



,-HfvL, . 

A < sfl ^££ji 



/-r—rHMp. , 

a < — ■^ L - 



(95) 



We now discuss the interaction contribution to the energy density. The expectation value of the 
energy-momentum tensor at leading order is given by 



(T^(x, r,)) = -2Re i / (T« (x, rj)H {I) 



(96) 



where is given in (|43[) . The energy-momentum tensor T$ is also in the interaction picture, and 
therefore, one finds that the energy density component is given by the Hamiltonian density in the 
interaction picture. Thus, we obtain 
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+ permutations . (97) 



The largest contribution in powers of j3r is of the form 
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where we have defined 

r), rf) = n(k hV ) - Q(ki,r/) . (99) 

The important point here is that at early times the oscillations due to the term ^^(^'V^vm) m equation 
(|98p severely damps this contribution, and at more recent times this already small energy density is 
further redshifted by the a(r])~ 4 factor. The same mechanism was also operating in the modified 
vacuum case discussed in [3]. 

We therefore conclude that the leading contribution from the interaction is in fact again given 
by the term proportional to 1/3^1 2 • By considering this term, changing the momentum integration 

variables as k{ — > yi = ^j^- and performing the relevant integral, we obtain the estimate 

\ / Planck Planck 

which does not strengthens the constraints coming from the free contribution. 

We also have to check the backreaction coming from the higher derivative interaction term. As 
before, the energy component of the energy-momentum tensor in the interaction picture is the Hamil- 
tonian density, in this case of equation (J75J). For subhorizon scales, the terms ip' and di<p scale in the 
same way, so the contribution to the energy density is 



(pW) « Re 



4A8a(7?) 6 J ( 2vr )l J ( 27 r)i jrfn i=1 

+ permutations . (101) 

By changing again variables as k i — > yi = kl7 ^ and performing the integration, we obtain the estimate 

(p( 3 ))«0/A 2 =eF 2 M p 2 ianck A 2 , (102) 

Observe however that already for A at the level of the supersymmetric GUT scale, A 4 > e-f^M 2 ^^, 
and therefore, once again, the constraints coming from the free contribution are not strengthened. 
The final result from equations ([32]) . (fSSjl . (fTUUD is 

M M 

W U \ < A P 2 ' anck ■ (103) 



6 Final results and conclusion 

Lorentzian symmetry could be broken at very high energies, for example by quantum gravity effects. 
The consequent modifications in the dispersion relations could leave detectable signatures in cosmolog- 
ical observables such as the temperature fluctuations of the CMBR, possibly allowing an experimental 
investigation of these aspects of the theory. 

In this article we have undergone a full general analysis of the effects of modified dispersion relations 
on the bispectrum, which is the leading contribution to the non-Gaussianities of the temperature 
fluctuations in the CMBR. We have in particular focused on dispersion relations that violate the 
adiabatic conditions at early times for a short period of time. The fact that the adiabatic condition 
is satisfied at the earliest times allowed us to fix unambiguously the initial conditions (vacuum). 
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The field equation in the presence of modified dispersion relations is difficult to solve, and therefore 
we have been using the WKB approximation for that scope, where possible. This approximation 
scheme is reliable and has been shown to be in quantitative, beside qualitative, agreement with the 
exact solutions, where available, see also [21]. It also allows to obtain general results. 

The universal features of the bispectrum for the modified dispersion relations have been discussed 
both in quantitative and qualitative terms. In particular, it has been shown: first, that when there 
is no WKB violation at early times, the result does not strongly differ from the standard slow-roll 
suppressed one. Second, that in the case of violation of the WKB condition at early times, the leading 
corrections to the bispectrum could be enhanced. The magnitude of the enhancement factors and the 
configurations for which these appear depend on the specific form of the dispersion relation. 

The largest enhancements for a given momenta configuration {k} arise if there exist a solution 
to the system of equations 

dySop({k},y)\ y=yt = S$ ({k},V*) + 0, 9?S of) ({k},V)\v=y. =0 Vra < n, (104) 
where the function S is defined in terms of the comoving frequencies oj as 



Sop(ki,k 2 ,k 3 ,y) 



( dy'r^u}(k h k m 1 ,y)-u(kj kJ-^U , h, j £ {1, 2, 3}, k in =max(k 1>2 ,3)- (105) 
3 Wj 7 



The schematic formula for the leading enhancements, taking care of the constraints (|103p from back- 
reaction, is 

• for cubic interaction as in the Hamiltonian (03]), see equations (|75p. 
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A\ ~ 1 r-r-;HM m . , x 

e]yu] rf^, (106) 



'0/3 

• for higher-derivative interactions as in the Hamiltonian (|78p . see equations (|87D . 

yl 




(107) 



We find differences when comparing our general results with those found by [3] using the modified 
vacuum approach to model transplanckian physics, which is based on imposing a cutoff on the momenta 
at a certain energy scale. The enhancements in the case of [3] arise only for the so called enfolded 
configurations, where the sum of the moduli of the external momenta vanishes, while for all other 
configurations the contribution to the bispectrum is the standard slow-roll suppressed one. Instead, 
in the case of modified dispersion relations violating WKB at early times, the enhancements could 
be enjoyed by different triangle configurations, depending on the particular dispersion relation at 
hand, leading to enhanced oscillations over different areas of the triangle space. In the lucky case, 
this could increase the detectability of the signal and possibly somehow alleviate the suppression due 
to the projection onto the two-dimensional surface following the decomposition of the bispectrum in 
spherical harmonics. We leave this point for future research. 

Also, the magnitude of the enhancements does vary with the form of dispersion relation. In the 
case of cubic coupling given by equation (fl3]) . obtaining the same magnitude as found in [3] would 
require quite particular conditions (that is, all the derivatives of the function in equation (|105|) to 
be zero on the enhanced configurations). On the other hand, for higher-derivative interactions as in 



3 Smaller enhancements are also possible for nearly critical boundary points, see footnote 1101 
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the Hamiltonian (|78p . the magnitude of the enhancements can be larger than the one found with the 
modified vacuum approach, if 

^) "y*>i- (108) 

For instance, for A ~ 10 3 H, e ~ \u\ ~ 1(T 2 , M™ . ~ 10 5 if, ^ ^ 0.5, n = 2 we obtain ~ 10 4 . 

3 " Ir I " Planck If* 3 i*V>n 

However, we stress that the presence of stationary points, in particular satisfying the condition 
(|108p . strictly depends on the form of the modified dispersion relation, and therefore could not be an 
easily occurring feature. Nonetheless, also in the worst case (n = 1 in (|107p : no stationary point), we 
can easily have enhancements because of boundary behaviour. For instance, in the condition of the 
example above they would be of the order 10 3 , which is of the same magnitude as those found in the 
modified vacuum framework for transplanckian physics. 
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A Appendices 

A.l Lorentz breaking action for the infiaton sector 

An action for the Lorentz breaking inflationary sector of our models can be written as follows: 

d 4 xy/^(C^ + C u ), (109) 



^ = -\g' LV d^d v <p - Y, c np D 2n 4>D 2 P4>, (110) 

n,p<n 

Cu = -A(0"%«„ + 1) - dxF^F^ = V^u v - V u U/1 , (111) 

where is the dynamical vector field necessary to properly define the Lorentz-breaking higher- 
derivative terms for the infiaton field (f>. Those terms are defined using the covariant derivative T>^ 
associated with what corresponds to the spatial metric seen by an observer comoving with Un, such 
that 

T>a(9v, + u li u v )=V a q lu , = 0, V 2n = (q^V u q^V a ) n . (112) 

In the action above, we have also constrained the vector field to have unit norm using the Lagrangian 
multiplier A, as it is usually done. 

If we choose a foliation of spacetime such that is orthogonal to the hypersurfaces of constant t, 
the modified dispersion relation obtained from the equation of motion following from the action (|109p 
is 

<Jp) = + 2£(-i) (n+p W£ 8 +p) • (H3) 

n,p 
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A. 2 A detailed example 



We consider now an explicit example of dispersion relation that violates the WKB condition at early 
times, given by 

' .2 _ P 4 „ ^ „. k 



(jj 



phys 



P" - x P < Pt 

C{p ~ Po) 2 + P>Pt' 



a{ri) 



This dispersion relation is identical to the Corley- Jacobson one with negative coefficient for p < pt [22] , 
but avoids the problem related with the presence of imaginary frequencies in the original proposal 
thanks to the modified behaviour for p > pt- We will however not discuss its phenomenological 
viability. 

The relation has a local maximum at p = p m = -A^ and is continuous and differentiable for 



P? 1 - 
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If we assume p(rj. 

p(v 3 



i %Pti the momenta p(j] 2 3 ) of horizon crossing (see figure [2]) are 
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(115) 
(116) 

(117) 



We also need to impose < ojq < \/2H which yields some conditions on c, e. We do not write the 
formulas as they are complicated and of little interest here. We make sure than in the following all 
the conditions are satisfied_L 



w p hys(p) 2 A 2 



A=10 3 H, c=499995.2 




Figure 2: Modified dispersion relation used for the example in section IA.2I 



We now solve the equation ([9]) as done in section 13.11 and compute the main contribution to the 
bispectrum for the cubic coupling in (|43p . coming from the integration over the interval [^n,^], as 
discussed in section 14.1.11 We neglect the other subdominant contributions. As expected, it is found 
that % (m ) ~ i] 2{3) . 

We compute the contribution to the bispectrum from the cubic interaction (|43|) using the formulas 
(1571) . ([59]) and ([63|) adapted to our case. We do not report their final expressions after being adapted, 

20 We require luq > in order to discuss an example with the good physical quality of not having a backreaction issue 
nowadays, which would place strong bounds on as we have discussed in section [5] 
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since the formulas are quite long and complicated but straightforward to be obtained. We find that in 
our example, the only enhanced configurations are close to the enfolded ones, where we have nearly 
critical points (see footnote fT0|) . We plot the bispectrum rescaled by the standard slow-roll result and 
the modulus of /3r in figure [3J taking A = lO^H. The dispersion relation itself is plotted in figure [2 

For comparison, we also plot the correction to the bispectrum for cubic coupling obtained in [3] 
using the modified vacuum approach to transplanckian physics (from equation (3.17) in [3]). We 
observe that the shape of the modulations and the magnitude of enhancement are different between 
that case and the one of the modified dispersion relation given by equation (|114p . 



10 3 H 
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Figure 3: Cubic interactions ~ £ 3 : a) Plot of the rescaled contribution to the bispectrum when all 
Green function are in the positive- frequency branch: negligible differences with the standard result 
(tiny superimposed oscillations), b) Plot of the rescaled contribution to the bispectrum when one of 
the Green functions is in the negative frequency branch: enhanced configurations for nearly critical 
points (see footnote [TUj) . c) Plot of the rescaled contribution to the bispectrum when one of the 
Green function is in the negative frequency branch in the case of the modified vacuum modelling of 
transplanckian physics [3J. 
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